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Abstract 

We make a proof -of -principle demonstration that improved energy resolution can be obtained in a negative-ion time 
projection chamber, by individually counting each electron produced by ionizing radiation. 
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1. Introduction 

Double beta decay is a rare nuclear decay mode, with 
a typical mean lifetime t > 10 19 y. It was recently ob- 
served in 136 Xe with a half-life t = 2 x 10 21 y [DEI- The 
zero neutrino mode of this decay, /5y8(Ov), appears the- 
oretically possible given that neutrinos have non-zero 
mass. Such a decay would violate lepton number con- 
servation and confirm the Majorana nature of the neu- 
trino. It has not been observed, despite significant on- 
going experimental effort world-wide. As in any rare- 
event search, substantial reduction of radioactive back- 
grounds are a high priority in the search for fj/3(0v). 
Given that the spectral expectation for /3/3(0v) is a sharp 
peak at the double beta decay Q-value, excellent en- 
ergy resolution can be used as an effective discriminant 
against those backgrounds which remain. But even in 
the absence of any radioactive background, the sensitiv- 
ity of a search for /3/3(Ov) may be limited by the high 
energy tail of the spectrum from the two neutrino de- 
cay mode, via the finite energy resolution of the detec- 
tor. The fraction of /3/3(2v) counts which contaminate 
the fi/3(0v) signal region is given by (3) 



where m e is the electron mass, Q = 2458 keV is the end- 
point energy of the decay, 5 is the FWHM energy resolu- 
tion at the endpoint and k is a coefficient which depends 
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weakly on 6. As an example, suppose a xenon detector 
obtains 6 = 0.05 (equivalently, a Gaussian <x = 2.1%). 
From [3], we find k - 1 . In a 20 tonne-year exposure 
(typical of a next generation experiment such as LZ20 
lED), 15 yS/3(2v) background events would be expected 
in the /3/3(0v) signal region. Larger values of 5 quickly 
lead to larger leakage. Only upper limits on the half- 
life for /3/3(0v) exist at present. This suggests that there 
is a significant motivation for developing a new type of 
xenon detector whose energy resolution approaches the 
statistical limit. 

In this work we make a proof-of-principle demonstra- 
tion that improved energy resolution can be obtained by 
individually counting each electron produced by ioniz- 
ing radiation. The idea has been previously suggested 
||5][6] but to our knowledge never implemented. Our de- 
tector is a negative ion time-projection chamber, shown 
schematically in Fig. [T] Rather than drift and detect 
the electrons created by ionizing radiation, we introduce 
oxygen as an electronegative dopant to capture the elec- 
trons. Oxygen ions 2 , rather than electrons, are then 
drifted to the anode. Just prior to arrival at the anode, the 
electrons are recovered from the O^ and amplified, us- 
ing a large electron multiplier (LEM) Q. Ideally, each 
electron is then individually counted. 

The dopant concentration and gas density are tuned 
so that electron capture occurs over a range of ~cm. 
We note that previous negative ion time-projection 
chambers have used either pure carbon disulfide (CS2) 
J8), mixtures of CS2 with noble gases [10) or ni- 
tromethane (CH3NO2) ifTTTl dopant. In this work, we 
employed a mixture of argon (66%, as a low-cost sur- 
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rogate for xenon), carbon dioxide (30%, as quench gas 
[ 12]) and oxygen (4%) as electronegative dopant. 

2. Apparatus 

The construction of our detector is shown schemati- 
cally in Fig. [T] An enlarged view of the region around 
the LEM is also shown, with a simulated electron 
avalanche of gain x30. The LEM was constructed 
from copper-clad G10 printed circuit board of thickness 
0.432 mm with 0.305 mm diameter holes. After drilling, 
a chemical etch was used to remove 0.023 mm of copper 
from the circumference of each hole, in order to increase 
the breakdown voltage. A segmented anode consisted of 
four concentric rings of radius 3, 6, 9 and 40 mm, with 
0.25 mm standoff between rings. Signals were digitized 
from the center three anodes, and the large outer ring 
was simply maintained at ground potential. The rings 
were used as veto to ensure complete collection of all 
electrons from each event. 

A charge-sensitive preamplifier with a gain of 
1.4 V pC l and recovery time constant t - 140 fis was 
connected to each ring, and the calibration was verified 
by introducing a 10 mV voltage step across a 1 pF ca- 
pacitance at the preamplifier input. Prior to introduc- 
ing O2 dopant, we first verified proper detector opera- 
tion and LEM gain under electron drift conditions in Ar 
(with 30% CO2 quench gas). The E = 5.9 keV gammas 
from the 55 Fe source produce no = E/w = 209 elec- 
trons, based on w — 28.2 eV. This value was obtained 
from a weighted average of the argon (w = 26.3 eV) and 
CO2 (w = 32.8 eV) constituent w-values fl4"l . a result 
which is confirmed experimentally [15|. These drift as 
a swarm to the anode, are amplified by the LEM and in 
the absence of O2 result in a single, fast output pulse 
from the preamplifier. Typical pulse rise times were 
measured to be O(10) ns. From the known quantities 
of initial and measured charge, the total effective LEM 
gain was determined. This is discussed further in Sec. 
[3] A flow-through system with proportional-integral- 
derivative feedback pressure regulation was used to mit- 
igate detrimental effects due to detector component out- 
gassing. Research-grade gas mixtures were used, and 
stable operation was observed on a timescale of sev- 
eral days. All results presented in this article were ob- 
tained at a gas pressure P - 0.250 ± 0.001 bar and 
T = 296 ± 1 K. 

There are three distinct regions of electric field in 
the detector: (i) the Zd = 110 mm drift region, with 
E d ~ 0.05 kV cm -1 , (ii) the z A = 0.432 mm elec- 
tron avalanche region E A in the LEM holes, with typ- 
ical maximum electric fields of several tens of kV/cm 



ifTBI and (iii) the Zi = 2 mm signal induction region, 
with Ei = 1 kV/cm. E A is a function of radial displace- 
ment from the center of a hole. Values quoted in this 
work are given simply by E A = AV/z A , where AV is the 
bias voltage across the LEM. The high electric field E A 
also provided the requisite energy to detach a captured 
electron from Or. The mechanism for this detachment 
process is thought to be collisional fPTl . 

The amount of oxygen dopant required to fully cap- 
ture all electrons created by 5.9 keV photons was first 
obtained from the measured reduction in pulse height, 
using electron drift operation. We found that 1 % O2 was 
sufficient to fully capture all electrons (though higher 
gas pressures required lower O2 concentrations). How- 
ever, the capture probability was sufficiently low that 
the resulting negative ion signal was distributed across 
about 10 ms. We found that with 4% O2, the result- 
ing signal was localized to within about 5 ms due to the 
more rapid capture of ionized electrons. This gave an 
average 22 fis between subsequent ion arrival times at 
the LEM, sufficient to allow most ions to generate a dis- 



tinct signal at the readout plane. In Sec. 4. 1 we quantify 
the fraction of missed ion pulses due to the proximity in 
time of successive pulses. In order to ensure negligible 
event pileup, the 55 Fe source rate was set at about 3 Hz, 
using a 0.88 mm diameter collimator. 

3. Data acquisition and analysis 

3.1. Data acquisition 

Negative ion drift event records of 25 ms each were 
acquired at 20 MHz, with 12-bit vertical resolution. A 
typical event obtained with an effective LEM gain of 
3 x 10 4 (AV = 880 across the LEM) is shown in Fig. [2] 
Note that effective LEM gain refers to the total observed 
electron gain from a single electron, and accounts for 
any reduction due to secondary electrons which did not 
reach the anode. The interesting portion of the event as 
digitized is shown in the lower panel. An enlarged view 
with the effect of the charge preamplifier decay time re- 
moved is shown in the upper panel. A total of n e = 157 
charge pulse steps were counted: 103 on the center an- 
ode, 54 on the adjacent ring anode and on the outer 
ring anode (not shown). Each step is interpreted as cor- 
responding to a single electron that was recovered from 
an 2 ion, and amplified by the LEM. The possibility 



of pileup is discussed in Sec. 4.1 



3.2. Analysis algorithm 

Our analysis algorithm employs a simple edge- 
detection scheme: for a given sample i, the average 
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Figure 1: (Left) Cross-section schematic of the detector, showing (1) initial generation of electron-ion pairs, (2) rapid electron capture on oxygen 
and (3) negative ion drift. (Center) Enlarged cross-section schematic of the LEM, showing a simulated electron avalanche (with gain x30) as 
obtained from the Garfield software 1 13]- (Right) Close-up photographs of the LEM (top) and segmented anode (bottom). 



value of the preceeding samples are subtracted from 
the average value of the following n a samples. The anal- 
ysis shown in this paper used n a - n\, - 10 samples. 
As might be expected, increasing the number of sam- 
ples contained in the average results in a higher sig- 
nal to noise ratio, allowing smaller pulses to be dis- 
cerned. We found as many as 5% more pulses with 
«o = nb — 50 samples. However, this comes at the 
expense of also increasing the number of pulses lost to 
pileup. In other words, the algorithm cannot resolve 
two single electron pulses that are separated in time by 
< HsI/adc = 1 yus. In this equation, n s = n a + tit, is 
the number of samples used to find the edge of a single 
electron pulse and /adc - 20 MHz. These pileup losses 
could be easily recovered by increasing /adc- 

The algorithm is capable of finding arbitrarily small 
pulses in an event record. In order to make a meaningful 
determination of the number of recovered electrons, a 
pulse height threshold must be chosen. A simple choice 
is the pulse height value at which the slope of the dis- 
tribution abruptly steepens, as indicated by the arrow in 
Fig. [3] (inset). This point corresponds to the onset of 
spurious noise pulses, and is determined by the input 
capacitance to the preamplifier. Accordingly, slightly 
higher thresholds were required on the larger-area an- 
odes. We note that Fig. [3] corresponds to the center 
anode. The thresholds employed in this analysis were 
[0.32 0.35 0.42] mV, for the center, ring and outer ring 
anodes. We also investigated a threshold determination 
based on the observed pulse height distribution in a re- 
gion where no signal is present, such as the initial ~ 1 
ms of an event record. Both methods lead to similar 
results. 

Additionally, we explored a slightly more sophisti- 



cated analysis algorithm in which the existence of a 
pulse in the event record was determined by the 2nd 
derivative of the raw event record. In other words, a 
pulse was identified by a zero-crossing in the double 
differential. As with the simple scheme, the pulse can- 
didate was still required to pass a threshold requirement. 
Similar results were obtained with both algorithms. 

We note that the analysis avoids double-counting an 
electron pulse which may have had its signal split across 
two adjacent rings. This scenario was found to occur on 
approximately 5% of pulses. 

3.3. Single-electron pulse height distribution 

Figure [3] shows the distribution of single electron 
pulse heights corresponding to the conditions exempli- 
fied in Fig. [2] The distribution of pulse heights resulting 
from avalanche multiplication is expected to be expo- 
nential 1 18 1 and is often referred to as a Furry distribu- 
tion. We observe a steeper rise at small pulse heights. 
The origin of this effect is a displacement of the z coor- 
dinate at which an electron avalanche was initiated. 

Consider the enlarged view of the LEM avalanche re- 
gion shown in Fig.[T] In that example, the avalanche ini- 
tiates at z ~ -0.15 mm. In the absence of (electron 
drift), the variation in the z coordinate of avalanche ini- 
tiation is cr ss 0.05 mm, as determined from our detach- 
ment model (described in the appendix). When 2 is 
introduced, avalanches are (on average) initiated farther 
into the LEM, because electrons must first be detached 
from 2 ions. This lowers the effective gain for a major- 
ity of electrons, as reflected in the reduced pulse height. 
This result is confirmed by our detachment model. A 
simulated pulse height spectrum obtained from the de- 



3 




0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3 

Time [ms] 



> 


50 


i i 


i i i 


i i i i i 

^enlarged view shown 


1 




i i i 




40 






1 above with amplifier 








CD 

-c 


30 






| decay (r) removed \ 








itu 


20 
















Amp] 


10 


















i i 


i i i 


1 1 1 1 1 


I 




i i i 



-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Time [ms] 



Figure 2: Example event record from 5.9 keV energy deposition in Ar-CC>2 with 4% O2 dopant at p = 0.25 bar. The LEM was operated with 
AV = 880. (Lower panel) The event record as digitized. (Upper panel) Enlarged view of the event record, after correcting for the charge 
preamplifier decay constant r = 140 /is. (Inset) Additional detail over a 60 /is window. 
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Figure 3: The measured pulse height distribution of single electrons 
on the center anode, as exemplified in Fig.[2](black steps). The func- 
tion y = exp(-x/g) was fit to the tail of the measured distribution 
(solid red curve, with small-x extrapolation shown dashed). The pulse 
height distribution obtained from the detachment model is also shown 
(dotted blue steps). Equivalent effective LEM gain is indicated along 
the top axis. (Inset) The pulse height distribution near pulse detection 
threshold (< 1 mV). 



tachment model is also shown in Fig. [3] normalized to 
the number of counts above the software threshold. 

3.4. Gain curves 

Gain curves for the LEM are shown in Fig. |4] for 0% 
(stars) and 4% (circles) O2. The addition of O2 lowers 
the effective gain at any particular voltage. With 0% O2 
the signal corresponding to the full-energy peak is a sin- 
gle pulse with a fast O(10) ns rise time. The gain curve 
for 0% O2 was obtained from that total pulse height, as 
described in Sec. [2] The curve for 4% O2 was obtained 
from fitting y = exp(-x/g) to the tail of the single elec- 
tron pulse height distribution, as shown in Fig. [3] The 
decay constant g of the exponential (in mV) was taken 
to indicate the mean gain at each voltage. Considering 
Fig. [3] it is expected that this procedure may result in a 
modest overstatement of the effective gain. We note that 
the present work does not require a precise knowledge 
of the gain; however, it is critical to verify that the gain 
is sufficiently large to see single electron pulses. 

4. Results 

A typical energy spectrum obtained from counting re- 
covered single electrons on the center and ring anodes 
is shown in Fig. [5] Software data selection cuts required 
< 1 electron pulse in the first and last 3 ms of events, 
and < 5 electron pulses on the outer ring anode. The 



energy resolution is discussed in Sec. 4.2 



Figure 4: Gain curves for the LEM. The dotted crossed lines indicate 
the data point corresponding to the examples shown in Fig. [2] Fig. [5] 
and Fig. [5] Where not visible, uncertainty is smaller than the data 
point. 



4.1. Counting losses 

Figure [6] shows the number n e of counted electrons as 
a function of LEM bias voltage. A maximum electron 
counting efficiency of 0.78 was obtained. We consider 
several likely electron loss mechanisms: 

(a) electrons which were not recovered from ■ 

(b) electrons which were recovered from Oj but did not 
produce an avalanche above threshold. 

(c) electrons which were recovered, produced an 
avalanche above threshold but were not resolved 
due to the finite value of /adc- 

(d) ions which collided with the top (in voltage) of 
the LEM, rather than being directed into one of the 
multiplication holes. 

(e) "gas-phase chemistry" |6|. An example would be a 
dissociation reaction which left an appreciable frac- 
tion of e.g. O , whose higher electron affinity then 
precludes any hope of recovering the captured elec- 
tron. Quantifying this effect is beyond the scope of 
this work, and is not discussed further. 

Addressing items (a) and (b), the trend of the data 
suggests that the full number of initial electrons could 
be recovered, if only AV could be increased by another 
~ 50 V. However, in this bias regime frequent electrical 
breakdown makes the LEM gain structure unusable. 

Addressing item (b) in particular, we can use our de- 
tachment model to make a rough estimate of the fraction 
of single electron pulses which were lost below the soft- 
ware threshold. Looking at the simulated curve in Fig. [3] 
(inset), the predicted fraction of pulse heights below the 
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Figure 5: The energy spectrum obtained from counting individual 
electrons recovered from Oj ions, obtained in Ar-C02 at p = 0.25 bar 
with 4% 2 dopant. The LEM was operated with AV = 880. The full- 
energy peak is 5.9 keV and the argon k-shell escape peak is 3.1 keV. 



Figure 6: Average number n e of counted electrons (with lcr width) as 
a function of bias voltage AV across the LEM. Corresponding Ea = 
AV/za is indicated along the top. 



software threshold is about 0.41. Under the experimen- 
tal conditions shown there (AV = 880), the total frac- 
tion of electrons not counted was 1 - 136/209 = 0.35. 
We see that the model slightly over-predicts the losses 
below software threshold. This is not surprising since 
there may be other factors (not addressed by the model) 
which influence the shape of the pulse height distribu- 
tion. Nevertheless, the rough agreement is encouraging. 

Addressing item (c), the fraction of pulses lost to 
pileup is approximately n c n s /N, where n c is the number 
of pulses occurring on a channel during the N samples 
spanning a typical event. As sho wn i n Fig. [2] typically 
N ~ 6x 10 4 samples. From Sec. 3.2 n s = 20. We must 
have n c < n e , where n e is the total number of pulses in 
an event. The limiting case occurs if all the electrons 
arrive at the same channel. Taking n c = 136 as in Fig. [5] 
this simple estimate implies that the typical maximum 
fraction of electrons lost to pileup is 0.04. The result 
is confirmed by a detailed analysis of the time between 
successive pulses in all events in our data sample. 

Addressing item (d), the next obvious parameter 
which might be increased is the induction field strength 
Ej. Previous studies ifTBI have shown that the trans- 
parency of a LEM for electrons depends on Ej. For 
LEM bias values AV < 880, we found that by dou- 
bling Ej from 1 kV/cm to 2 kV/cm we obtained a 20% 
increase in the number of recovered electrons. Fur- 
ther increase of E, did not yield additional recovered 
electrons. Interestingly, when operating at higher bias 
values AV > 880 with Ej = 2 kV/cm, the maximum 
number of recovered electrons was still limited to about 
xO.78. 



An additional complication of operation at E, > 
1 kV/cm is a distortion of the pulse shape. With Ej = 
1 kV/cm individual pulse shapes were as expected: a 
voltage step with fast O(10) ns rise time, followed by 
recovery with r = 140 With Ej = 2 kV/cm, we 
observed the same initial fast rise time, followed imme- 
diately by a slow ~ 3 /-is rise time, followed by recovery 
with the expected time constant. The step height of the 
slow rise time was less than half the step height of the 
fast rise time, and was observed even in the absence of 
oxygen. This leads to the interpretation that it is due to 
positive ion movement between the lower (in voltage) 
electrode of the LEM and the anode. Physically, this 
implies that the electron avalanche continued for sev- 
eral hundred microns beyond the lower electrode plane 
of the LEM. 

We have assumed that w-value is unaffected by the 
presence of the oxygen dopant. This seems reasonable 
since oxygen serves to remove, rather than add, ionized 
electrons. However, for completeness we note that if the 
w- value were higher (lower) by e.g. 2 eV, the predicted 
resolution shown in Fig. [7] would be lower (higher) by 
only a few percent. This is a small effect - roughly 
equivalent to the line width of the curve. At the same 
time, the expected maximum n e shown in Fig. [6] would 
decrease (increase) by about 7%. 

4.2. Energy resolution 

The best energy resolution we achieved under elec- 
tron drift conditions (in the absence of O2) was cr/p = 
0.11. In Fig. [7] we show the energy resolution o~/n e as 
a function of n e , as obtained from negative ion drift. 
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Statistical uncertainty is indicated by the thin verti- 
cal bar, and in several cases is smaller than the data 
point. Sources of systematic uncertainty, indicated by 
the larger vertical rectangle, include (1) varying n a j, be- 
tween 10 and 20 samples, (2) variation due to binning 
and fitting and (3) differences in results between two in- 
dependent analyses. The last item implicitly accounts 
for variation in the determination of the pulse detection 
software threshold. 
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Figure 7: Observed energy resolution (circles) compared with the ex- 
pectation for a binomial process (thin solid), the prediction from Ref. 
1 6 1 and our prediction as given by Eq. [2](thick solid). Statistical uncer- 
tainty is indicated by the thin vertical bars, and systematic uncertainty 
(discussed in the text) is indicated by the vertical rectangles. 



The expected detector resolution for electron count- 
ing is 



<r/n e 



F | (1-p) 
n n p 



(2) 



where F = 0.20 |[T9l is the Fano factor in gaseous ar- 
gon. This equation states that the resolution depends 
on the Fano fluctuation in initial number no of ionized 
electrons, and on the binomial process by which elec- 
trons are recovered from Oj and amplified. As men- 
tioned previously, we define n e = nop, where p is the 
total probability for an electron to be counted. 

The agreement between Eq. |2]and the data is not sat- 
isfactory for any value of n e . We are therefore lead to 
consider other processes which remain unaccounted in 
our analysis. We note that if the measured n e values 
were ~ 30% smaller, our data would appear to agree 
well with Eq. [2] This leads to the suspicion of a spuri- 
ous double-counting of recovered electrons. We have al- 
ready removed the possibility of double-counting pulses 
across adjacent anodes. There remains the possibility 



of double-counting due to some sort of feedback mech- 
anism. For example, we know that copious UV pho- 
tons are generated during the electron avalanche. Could 
there be a small probability that one of these photons 
eject an electron from e.g. the wall of the LEM, lead- 
ing to a spurious electron pulse? We have investigated 
this possibility and found no concrete evidence to sup- 
port the hypothesis, but this does not allow us to rule it 
out. For completeness we note that the discrepancy we 
observe cannot be accounted for by any reasonable vari- 
ation in the Fano factor. We also note that our data do 
not support the theoretical energy resolution proposed 
in El. 

5. Summary and future work 

We have demonstrated that improved energy reso- 
lution can be achieved in a negative ion time projec- 
tion chamber by individually counting the electrons pro- 
duced by ionizing radiation. The initial results are en- 
couraging, but many questions remain to be addressed. 
Higher electron gain and a higher electric field Ea are 
clear priorities. Future work will study the electron 
counting performance using proportional wire and MI- 
CROMEGAS gain structures, as well as optimal gas 
mixtures (including our desired target isotope, xenon) 
and operation at increased gas pressure. 

Appendix: a simplified electron detachment model 

In order to further understand the physics underlying 
the single electron pulse height distribution discussed in 
Sec. 3.3 we employed a simplified model of the de- 



tachment process. It was assumed that ions drift in 
an ideal gas of neutral argon molecules in thermal equi- 
librium. The ions make elastic collisions described by 
a cross section calculated assuming a long-range polar- 
ization interaction with the neutral atom 1 20 1 . They are 
accelerated by a constant electric field between colli- 
sions. The electron detaches from the ion when the cen- 
ter of mass energy of the ion-neutral system (between 
collisions) exceeds a threshold s D . Using momentum 
transfer theory [2111 . the model produced a lifetime for 
electron detachment for a given electric field strength 
and so, and a distribution of detachment probabilities 
for a given electric field profile. It would seem natu- 
ral to associate e d with the electron affinity, which is 
e = 0.45 eV l22l for O2. Based on experimental deter- 
minations of the collision cross section for detachment 
of both ions [ 23 1 and several types of halogen ions 
11241 . it appears that can often exceed the electron 
affinity by x2 or more. 
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Figure 8: Simulated distribution of z coordinate at which electron 
avalanches begin, with (solid) and without (dashed) Ot. The verti- 
cal dotted lines indicate the physical extent of the LEM. 

We therefore treated ed as a free parameter. Different 
values of e D resulted in different z detachment profiles. 
These were incorporated into Monte Carlo simulation 
to study the effects of ion detachment on single-electron 
pulse height distributions. The simulated pulse height 
distribution for the avalanche gain for eq — 3 eV is 
shown in Fig. [3] The point of avalanche initiation which 
resulted from the detachment profile (also obtained with 
Ed — 3 eV) is shown in Fig. [8] The results were obtained 
for a gas mixture of ArCC<2, using Garfield++ ll25ll in 
conjunction with Gmsh 11261 and Elmer [27]. The sim- 
ulated distribution for electron drift is also shown. We 
note that by default, the predicted gain was about xlO 
too low. This may be partially due to Penning processes 
which are not accounted for in the Garfield++ simula- 
tion ll28l . As an interim solution we explored the ef- 
fect of varying the default Townsend coefficient data in- 
put to Garfield++, and found reasonable agreement (as 
shown here) with a uniform increase of 19%. This is 
well within the range allowed by experimental uncer- 
tainty 12911301 . 
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